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In thin films of smectic-C liquid crystals, localized regions containing additional smectic layers 
form circular inclusions that carry a topological charge. Such inclusions nucleate a companion 
topological defect. These inclusion-defect pairs are modeled as topological dipoles within the context 
of a one-coupling constant approximation to the 2D Frank free energy. Deviations of the dipole 
direction from a preferred orientation cause the dipoles to acquire a logarithmic charge. Thermal 
fluctuations of the dipole direction are calculated and found to be large, scaling as the logarithm of 
the system size. In addition to dipole-dipole interactions arising from the topological charges, we 
also find that the thermal fiuctuations of the dipole directions are coupled through a preference for 
global charge neutrality of the logarithmic charges. 
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I. INTRODUCTION 

Phases of matter with broken continuous symmetries 
exhibit topological defects, which can be thermally ex- 
cited or induced via external stresses. The surface inter- 
actions of a foreign inclusion within an ordered medium 
can induce nontrivial distortions in the bulk medium and 
lead to the nucleation of defect structures Even in 

cases where no defects are present, the distortions that 
arise can often be described by virtual defects within the 
foreign inclusion itself in much the same way that elec- 
tric fields outside a conductor can be described by image 
charges within the conductor. Two such experimental 
systems which have recently been studied are Langmuir 
monolayers and nematic emulsions [0J|]. 

At sufficiently high density, a Langmuir monolayer 
forms an orientationally ordered liquid-condensed phase 
in which the order is characterized by a 2D unit-vector 
order parameter, which we call the c-director Un- 
der appropriate conditions a liquid-condensed region can 
contain isolated domains of the rotationally isotropic 
liquid-expanded phase, that is, regions wherein the ori- 
entational order is absent. Along the boundary separat- 
ing the two phases, there is typically a preferred orienta- 
tion of c relative to the tangent vector to the boundary. 
The inclusions (liquid-expanded regions) can thus lead to 
a distortion of the order-parameter field of the ordered 
medium (liquid-condensed region) surrounding them. 

Nematic emulsions are dispersions of droplets of ne- 
matic liquid crystal in an isotropic host Inverse ne- 
matic emulsions consist of surfactant coated spherical 
water droplets dispersed in a nematic host The 
surfactant can be prepared such that the nematic direc- 
tor prefers to align either parallel or perpendicular to a 
surfactant coated surface. Thus the droplets distort the 
local nematic field and, in the case of normal bound- 
ary conditions, the droplets themselves carry a topolog- 
ical hedgehog charge which must be compensated for by 
creating additional defect structures within the ordered 
medium. 

Experimental observations (l||^J^ in both Langmuir 
monolayers and inverse nematic emulsions have been ex- 
plained theoretically with some success ^,|[. Neither, 
however, is the simplest system one could imagine for 
the study of the effects of foreign inclusions in an ordered 
medium. In Langmuir films the line tension of the inter- 
face between the two phases is so low that orientational 
order in the surrounding medium and the shape of the 
inclusion are strongly coupled, greatly complicating the 
theoretical treatment. In nematic emulsions, the elastic 
theory describing director distortions is highly non-linear, 
and as a result, it is very difficult to find exact solutions 
for the director around a droplet-defect pair. In this pa- 
per, we report observations of freely suspended inclusions 
in another ordered system, smcctic-C films, and we de- 
velop a theory for their energetics. 

There are several types of circular inclusions with fixed 



boundary conditions that can be inserted into a freely 
suspended smectic-C film. The type we are concerned 
with here are islands, places where the film contains an 
integral number of additional smectic layers (Fig [l]) . Is- 
lands can be introduced to a freely suspended film by 
(1) rapidly reducing the film area, (2) blowing across the 
film to pull material from the meniscus region or (3) heat- 
ing the film near the smectic-to-isotropic or smectic-to- 
nematic phase transition, at which point isotropic or ne- 
matic droplets form in the film and remain as islands af- 
ter cooling. Typical islands have radii from 50 to several 
hundred microns. The boundaries of these inclusions are 
smectic dislocations, and the c director has a preferred 
orientation relative to the tangent line to such disloca- 
tions fioll thus providing our fixed boundary condition. 
This fixed boundary condition allows us to consider the 
interior and exterior regions of an island as separate sys- 
tems, analogous to the way a conducting surface allows 
one to independently consider the electrostatic config- 
urations interior and exterior to the conducting surface. 
While island-defect pairs are relatively easy to create and 
image (Fig. |l|) , they have the disadvantage that often the 
islands shrink, and it is difficult to control their size. This 
makes the creation of several islands of uniform size very 
difficult. 

A second type of inclusion in smectic-C films are 
isotropic or nematic droplets, found in materials that 
have a direct phase transition from the smectic-C to the 
isotropic phase. Droplets have the disadvantage com- 
pared to islands that as one nears the isotropic transition, 
the birefringence of the film diminishes, and it becomes 
difficult to image the c-director. A third type of inclu- 
sion consists of the direct insertion of small (10 micron) 
glass spheres into the film. While this in principle is the 
best technique, since size can be accurately controlled 
and the birefringence is unaffected, in practice it is dif- 
ficult to maintain glass spheres larger than 5 microns in 
the film. Larger spheres either rupture the film or simply 
pass through. 

Regardless of its type, a circular inclusion with rigid 
boundary conditions along its perimeter carries a topo- 
logical charge of +1. Experimentally it is observed that 
each island nucleates a companion —1 defect out of the 
bulk, which maintains a preferred separation from the 
inclusion while undergoing visible radial and azimuthal 
orientational fiuctuations with respect to the center of 
the inclusion. 

In this paper, we present a theoretical study of these 
circular inclusions in smectic-C films. We also report 
briefly on preliminary experimental observations that are 
in qualitative agreement with our theoretical analysis. 
The theoretical treatment of circular inclusions with rigid 
boundary conditions in smectic-C films is simpler than 
the treatment of similar defects in Langmuir films and 
nematic emulsions since shape changes in the bound- 
ary present in Langmuir films do not occur and within 
the one coupling constant model of smectic-C elastic- 
ity we avoid the nonlinearities inherent in the treatment 
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of 3D nematics. In fact, the equation describing the 
order is simply Laplace's equation, = 0, where 

c = (cos$,sin$). Because we impose rigid boundary 
conditions, our 2D results are very similar to those ob- 
tained in SD nematic emulsions with similar boundary 
conditions. Fluctuations, however, are more pronounced 
in 2D than in 3-D. 

If the angle $ is constrained to make a constant angle 
(j)a with respect to the tangent to the boundary of an 
inclusion, then an inclusion necessarily creates (or "is") 
a topological defect - a vortex, or equivalently, a discli- 
nation with positive topological charge. If the bound- 
ary conditions at infinity enforce alignment along a spe- 
cific direction the total topological charge must be zero, 
and, as a consequence, a companion defect with negative 
charge must be created for every inclusion. The inclusion 
and its companion defect form a topological dipole that 
interacts with other distant inclusion-defect pairs via a 
2D dipolc-dipolc interaction. 

In 31? nematic emulsions, fluctuations in both the dis- 
tance between a droplet and its companion defect and 
the angle 6 between the droplet-defect dipole and the 
direction of the far-field alignment are negligible, with 

where T is the temperature, X is a Frank elastic constant 
and a is the droplet radius (~ 1/im). Since K « ksT/l, 
where I is approximately a molecular length ( 300nm) we 
find that (((56')2) « l/a « 0.003 « 1. 

In our 2D analog, the Frank elastic constant K2D has 
units of energy and is typically of order dK^u, where 
d is the film thickness. Taking d ~ nl, where n is the 
number of layers in the film, we then find {{5rd/rdf') ~ 
kBT/K2D ~ 1/n. Our calculation actually leads to the 
numerically smaller result {{Srd/rd) ) ~ 1/8?™. In ad- 
dition, logarithmic singularities in two-dimensions lead 
to a logarithmic dependence on the system size, R, in 
({Sq^y) « log (R/a) {kBT/2nK2D)- ^ 

Laplace's equation in two dimensions admits two fun- 
damental types of point defects: a topological defect 
with (f> — t'An^^{y/x) and a logarithmic defect with 
(j> = log(a;'^ + y'^). In what follows, we will refer to the 
charge of the topological defect as the topological charge 
or magnetic charge and that of the logarithmic defect 
as the electric charge since the angle field of the latter 
is equivalent to the potential produced by a point elec- 
tric charge (see Appendix). To describe arbitrary ori- 
entations of the topological dipole with respect to the 
c-director at infinity, both kinds of defects are needed. 
The "charges" of the topological defects are fixed to be 
integers whereas those of the logarithmic defects are pro- 
portional to the angular deviation of the dipole from its 
favored direction. 

As in 3D nematic emulsions Q], widely separated 
inclusion-induced topological dipoles interact via a 
dipole-dipole potential arising from the arrangement of 



the magnetic charges associated with each inclusion. 
When the dipoles are not properly aligned with respect 
to the c-director at infinity, the inclusions will also carry 
a certain amount of electric charge, leading to an ad- 
ditional interaction. In three dimensions, the analogous 
interaction is unimportant because of the prohibitive cost 
of deviations from the direction of preferred alignment. 
In two dimensions, deviations can be large: this interac- 
tion is important, and it can be calculated. It predicts 
that correlated rotations of dipoles in which the total 
electric charge is zero are energetically less costly than 
those with non-zero electric charge. It should be possible 
to test this prediction experimentally. 

Our preliminary experimental observations are in qual- 
itative agreement with these predictions. Associated 
with each island is a companion negative defect (Fig. ||). 
The position of this defect fiuctuates visibly with angu- 
lar fiuctuations apparently greater than the radial fluc- 
tuations. Furthermore, defect fluctuations appear to be 
correlated when more than one island is present. Further 
study will be needed to provide unambiguous confirma- 
tion of these observations. 



II. THE MODEL AND SOLUTION 

Our free standing smectic-C film is closely approxi- 
mated by a 2D XF- model A typical film has a 
small number of layers, N, (roughly 4-10) which can be 
determined exactly with laser reflectivity measurements. 
The projection of the nematic director into the plane of 
the film provides us with our unit-vector order parame- 
ter c (the tilt is relatively constant over the area of the 
film). The effective free energy of these films contains 
two elastic constants - one for splay distortions (Ks) and 
another for bend (Kb)- 

P= \ I [^s(V • c)2 +Kb{V -K cf] d^x. (2) 

However, when the difference between Kb and Ks is 
small, the 2D XY-vaodel provides an accurate approxi- 
mation, 

/ [(V-c)2 + (Vxc)2]d2a; (3) 

= f / {y^?<fx, (4) 

where it is understood that the integral is taken over 
the region where the order parameter c exists and is 
nonzero. 

Let us now examine an ideal case of a circular inclusion 
immersed in an ordered medium described by an XY 
order parameter c = (cos sin <i>). We take Eq. (||) to be 
our free energy, where K is the Frank elastic constant and 
the integral is to be taken over the region that the ordered 
medium occupies (external to the inclusion and internal 
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to some boundary farther away), namely ft = Br — Ba, 
where B^ is a disc of radius r centered about the origin, a 
is the radius of the inclusion and R is the system size. We 
do not include any surface energies since we will assume 
the properties of the surface will remain fixed (i.e., strong 
pinning boundary conditions keeping c fixed and a very 
large line tension keeping the shape fixed). 

We will assume that the far field is uniform at the outer 
boundary taking $|r=_R = 0. Next we take the boundary 
condition at the surface to be of the form ^\r=a = 0+ f + 
4>a, where (j) is the azimuthal angle relative to the center 
of the disc and (/jq is the angle that the director c makes 
with the tangent vector to the boundary. This condition 
gives the inclusion a topological charge (winding number) 
of +1. The far-field boundary condition requires that 
the total charge be 0, thus requiring the existence of an 
additional topological defect with charge —1. We will 
find that this defect will sit fairly close to the inclusion 
and that the pair will constitute a dipolc, which can in 
turn interact with other such pairs via a dipole-dipole 
interaction. We will note that the value of (j>a will only be 
relevant if one knows the actual direction of the far field, 
and then it will determine the relative angle between the 
far-field direction and the direction of the dipole. 

Taking the radius of the defect to be ^ <C a and as- 
signing a core energy to the defect Ccorc ~ t^K we 
simply need to find solutions for $ that depend on the 
position of the defect and to extract an effective free 
energy, F[rd], as a function of the defect position r^. 
Minimization of Eq. (^) leads to Laplace's equation, 
= throughout except at points occupied by 
defects. The angle field of a single charge defect lo- 
cated at {xo,yo) = (rocos0o,''osin0o) = {ro,4>o) = is 
simply 



<^ro(r) = (l>roix,y) = Im(log|z - zq\) , 



(5) 



where z = x + iy. We note that (pro (r) satisfies Laplace's 
equation everywhere except at Tq. The angle field of a 
charge-g defect at Tq is q(f>ro- In particular, a charge —1 
defect has a field —(pro- 

There is an image solution to our problem of a disc 
with strong pinning boundary conditions and a uniform 
director field at infinity |13|,p[: 



i>Q +TT + 2k7T, 



(6) 



where r™ = (a^/rp, (f>o) and k is an integer, which may 
vary depending upon the point on the circle r = a at 
which we are evaluating the function. The function 



where 



~ 2 ~ 



(7) 



(8) 



has the following properties: $|r=a = {(t>+'^+4>a) niod 2n 
and $ — > as r ^ oo. Thus we can satisfy boundary con- 
ditions at r = a and at oo by creating a real —1 defect 



at (rd,0j]) in the medium outside the circular inclusion 
and a +2 defect at r = and a —1 defect at {a'^/rd, 0^) 
inside the disc. Since in the absence of boundary condi- 
tions at R, there would be a single -1-1 defect at r = 0, 
we may regard one -1-1 defect at r = and the interior 
— 1 defect as images. Even though this solution does not 
strictly satisfy our boundary condition at R, we claim 
that it is good enough for R significantly larger than a. 
To satisfy both boundary conditions would require an in- 
finite number of images inside a and outside R, but the 
additional images would have little effect on the energy 
as the interested reader may easily confirm. 

The collection of three charges can be viewed as a topo- 
logical dipole with moment p pointing from the induced 
defect in the medium to the center of the disc: 



-p{cos(j)d,sm(t)°d), p = rd + 



(9) 



The dipole makes an angle of tt -I- 0° with the direction 
of c in the far field. Sample configurations are shown in 
Fig (|^) . If c is parallel to the x-axis at infinity, then the 
induced defect lies above the disc, and hence p points in 
the —y direction when (j)a = Q (tangential boundary con- 
ditions); and the induced defect hes to the left of the disc, 
and hence p points in the -\-x direction when 4>a — —■^ 
(normal boundary conditions) . An intermediate configu- 
ration with 0a = is also shown in Fig. (||). 

The solution we have constructed is valid for any Vd^ 
and we now calculate the free energy of a single dipole, 
Fi[(a,rd,4>^d\^ ^ ^ function of rd at the fixed orientation 



Fi[{a,rd,(tPd)] = -TTK\og{ard)^2TTK\ogrd 



(10) 



—ttK log ^ — TiK log [rd + 2itK log + 

V TdJ 

A «^ a* \ , a 

= -ttK log — J ] + ttK log 7 + Cc 



We see that the preferred position for the defect is rj] = 
V2a. Expanding about this minimum with Srd ^ rd-^ r'^ 
we find. 



F[{Srd, 03)] ~ const + AttK(^ 



Thermal fluctuations in rd then satisfy. 




(11) 



(12) 



in qualitative agreement with our estimate from sec- 
tion (I) . Noting that this K is the two-dimensional 
constant K2D from section (I), we recall the estimate 
ksT/K « 1/n where n is the number of layers in the 
film (typically 4-10). This would suggest that the radial 
fiuctuations should be of order r^/lO for a 4 layer film. 
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Our experimental observations are slightly larger than 
this. 

We would now like to relax our assumption about the 
azimuthal location of the defect and calculate F[a, r^] 
where = (r^, (fid) and (j)d = 4>d + Hd- Equivalently, we 
now let the angle of the dipole (p = — p(cos 0^, sin 
vary as well as its magnitude (Fig. In doing so, we will 
in fact see that 0j] is the equilibrium azimuthal location 
of the defect. To find a solution, we have only to add 
to <I> in Eq. (|^) a solution to V^<i> = that is zero at 
r = R and equal to Sipd at r ~ a. Such a solution can be 
constructed from the electric-charge field Re log z — log r, 
which is constant for constant r. We find 



where 



$ = 20-(/.r, -0r--(7log(^-j, 



log(f) 



(13) 



(14) 



is the electric charge. The energy associated with this 
field is 



F[ia,rd,cl)d)]^F,[ia,rd,cf°d)] + 



ttK 
1^ 



bdf. (15) 



From this we can calculate the magnitude of the orienta- 
tional fluctuations: 



md 



2nK 



log 



R 



(16) 



Thus, in an infinite sample {(54>dY) is infinite, and there 
is no well-defined direction of p relative to the direction 
of c at infinity. In laboratory samples one typically finds 
R/a « 10^ - 10^ giving log(i?/a) « 5 - 9 and hence 
fluctuations in the angle are reasonably enhanced by this 
logarithmic factor. Finally again using our estimate of 
ksT/K l/n with n = 4 we flnd that for R/a = 10^ 
we have {{6(pd)'^) ~ 0.25. This corresponds to angular 
fluctuations of approximately 30 degrees, which should 
be far more noticeable than the radial fluctuations Srd- 
Observed angular fluctuations indeed appear to be of this 
order. 



III. MULTIPLE DOMAINS 

A single inclusion distorts the c-director field as we 
have just shown. If there are two or more inclusions, the 
distortions they induce will give rise to effective inter- 
actions among the inclusions. If the separation between 
inclusions is large compared to their radii, each inclu- 
sion behaves like a dipole, and there will be a pair- wise 
dipole-dipole interaction between them. In addition, if 
the dipoles are not properly oriented with respect to the 
far-field c-director they spawn "electric" charges within 



the inclusion, and these give rise to a 2D Coulomb inter- 
action between the inclusions. More importantly, since 
the values of these electric charges are not fixed, the usu- 
ally constant "self-energy" contribution will not be con- 
stant and will greatly favor total charge neutrality within 
the system. This in turn will imply that fluctuations of 
the dipoles adhering to electrostatic charge neutrality will 
be enhanced over those that do not. 

Before proceeding we wish to clarify our notation for 
a single inclusion-defect pair. For simplicity and with- 
out loss of generality, we will continue to assume that 
the far-field c-director points in the e^, direction. From 
the previous section we conclude that a single inclusion- 
defect pair can be fully described by [r,p, a, (j)a], where 
r is the position of the center of the inclusion (no longer 
required to be at the origin) , a is the radius of the inclu- 
sion, 0a tells us what the boundary condition is at the 
inclusion boundary, and 



-p(cos 0d,sin (fid); 



bd = 4'°d 



TT 
2 



P = rd-\ ; rd = r°^+Srd. 

Td 



(17) 
(18) 

(19) 



via q 



We then found it convenient to code 
5(t)d/\og[R/a). 

Now we would like to calculate the free energy for a 
system of such inclusion-defect pairs, F([ri, p^, a^, (/)ai]), 
where we assume that this description of multiple inclu- 
sions is valid when R^ jr^ — r^l a for all pairs (i, j). 
As we are now aware, each isolated inclusion-defect pair 
can be thought of purely in terms of the external topo- 
logical defect and an arrangement of topological and log- 
arithmic defects inside the inclusion. Of course, when 
there is more than one inclusion, additional images must 
be introduced in the interior of each inclusion to sat- 
isfy boundary conditions. The full calculation of the en- 
ergy is then reduced to calculating the energy from all of 
the topological charges and electric charges. As long as 
the inclusion-defect pairs are reasonably separated, how- 
ever, one can obtain a good estimate of the energy by 
ignoring the contributions made by the additional image 
charges. The energy decomposes in a natural way into 
three pieces: the first is the contribution each individ- 
ual inclusion-defect pair makes when its electric charge 

= 0, the second is the contribution from the dipole- 
dipole interaction among the topological defects, and fi- 
nally, the third is the total electrostatic energy arising 
from the electric charges: 

F({[rj,pi,ai,0ai]}) ~ '^Fi{ai,rdiAai) 

i 

+ 2ttK ( '^^ - 2(p, ■r.jQJpj ■r^j) \ 

i<j \ ^ij / 



6 



- ^ TiKqiqj log 



(20) 



i<3 



Yj. The 



where Fi is defined by Eq. ( pl)|) and Vij = 
second term is the dipole-dipole energy and the third and 
fourth are the electrostatic energy. Again we emphasize 
that this is only an approximate form for the energy valid 
when the pairs are sufficiently for apart. 

To see the importance of the first electrostatic 
term in the energy we now turn to a more spe- 
cific example. For simplicity we will assume that 
our system consists of only two inclusion-defect pairs: 
[0, pi, a, — 7r/2], [r^e^;, P2, a, —tt/2]. We have taken the in- 
clusions to have the same radii and to have the same 
normal boundary conditions. In addition we have placed 
one inclusion at the origin and the second at r^, where 
Ts is parallel to the preferred orientation of the inclusion- 
defect dipoles (Fig. (|^)). Assuming that the magnitudes 
of the dipoles are fixed at the value for a single inclusion, 
rj], we can write 



p(cos 0i,sin 



(21) 



where p = a{\/2 + 1/V2) and = n + (pd (note that 
(j)'^^ = —IT so that (fif ~ Q). This implies each inclusion- 
defect pair carries an electric charge 



Qi = 



Using Eq. (E^) we can calculate 



F[qi,q2,rs] = 2Fi[(a,rrf,- 



(22) 



(23) 



2TrK 



P1-P2 2(pi • rs)(p2 • r^) 



nK{qi + q^f log (^] ~ li^Kq^q^ log (^) 



where again the first term is simply twice the energy 
from Eq.([lo|) for a single inclusion with ^(i)^ = 0. Taking 
9+ = (ft + 92)/2 and q_ = (qi - q2)l2 = qi - q+, we can 
simplify now to find 



F[q+,q_,r,]^2F,[ia,rd,' 



(24) 



p f f R 

-2t:K— cos 2q+ log — 
ri \ \ a 



2^Kql{\og{^\+\og(^ 



2t: Kql log (^) 



We see immediately that qj^ — q^ = is the minimum, 
which corresponds to c/ji — (f>2 — where the dipoles are 
aligned with the -l-x-axis. More interestingly we see that 
g+ and g_ are decoupled. The equipartition theorem 
then implies 



ill) 



M(fJ(log(f))%log(f)+log(, 



{Q-) = 



inK log {^) 
For the (i, 's we have, 



(25) 



(26) 



{S^D = {S4) - {ill) + id)) (log (f ))' (27) 

{S<t>M = {{ql) - (d)) (log (I) ) ' . (28) 

In the limit R — > cx), fiuctuations in q^ tend to zero. The 
ratio of 6(l)i to qi, however, diverges so that in this limit. 



mi + 6hf) 



AttK 



is finite and 



keT \og{R/a) 
AttK log(rs/a) 



(29) 



(30) 



diverges. Note that fluctuations in (((50i — &412Y) in the 
two inclusion system are much larger than fiuctuations 
in 5(j)d for a single inclusion. Of course, the ratio R/ a in 
a real system is such that log(i?/a) is not much greater 
than 1. Let us compare {54'1) and (i50i(502) to {{S(f>ci)'^) 
from equation (|6|) for realistic values R/a = 10'^ and 
rs/a = 10. We find that {S(f>^) « {1.75)md)'^) and 
that {6(j)i5(t)2) w -(1.25)(((50d)^). For n = 4 this would 
correspond to angular fluctuations of a single dipole of 
approximately 40 degrees versus angular fiuctuations of 
only 30 degrees for an isolated domain. Furthermore we 
see that the fiuctuations are highly correlated. Thus we 
see that the presence of multiple domains can enhance 
the magnitude of the fiuctuations a noticeable degree. 

Had we not taken the separation vector to be in the di- 
rection of the far-field director we would have first been 
obliged to find the new equilibrium values for and 
g_. Not surprisingly we would have found that the g+ 
equilibrium value would no longer be zero, however the 
term in the energy for q- would remain the same. Since 
the fiuctuations in g_ dominate the fiuctuations in g+ 
we would still expect the above observations concerning 
the magnitude of the fiuctuations to be fairly accurate. 
Furthermore, in the presence of many domains, the corre- 
lations between them may be quite complicated but we 
will still expect some enhancement in the angular fiuc- 
tuations arising from those fiuctuations which preserve 
electric charge neutrality. 



IV. SUMMARY AND CONCLUSIONS 

Smectic-C films are adequately described by the 2D- 
XY model. Circular inclusions in the film with appro- 
priate boundary conditions carry a topological charge. 
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Global boundary conditions require the inclusions to nu- 
cleate companion topological defects out of the surround- 
ing ordered medium, providing a good method for study- 
ing the properties of a small collection of defects. Herein 
we have described how director configurations can be de- 
scribed systematically by externally nucleated topologi- 
cal defects and associated image defects lying within the 
inclusions, provided we allow for logarithmic as well as 
topological defects. The electric charge required is deter- 
mined by the orientation of the topological dipole created 
by the inclusion and its companion and by the direc- 
tion of the far-field director. The value of the electric 
charge fluctuates producing visible fluctuations of the 
position of the companion defect about its host inclu- 
sion. When multiple inclusions are present they inter- 
act through their topological and electric charge distri- 
butions. The topological dipole distributions associated 
with each inclusion-defect pair lead to a dipole-dipole in- 
teraction between them an to an alignment of the dipoles. 
However, thermal fluctuations of the dipole direction are 
large (the fluctuations the companion about the host) 
and carry electric charges with them. Orientational fluc- 
tuations of the dipoles are coupled via the interaction of 
these electric charges. By far the strongest correlations 
arise from the overriding desire for global electric charge 
neutrality in the system (the energy of an electric charge 
in 2D diverges with the system size). We have presented 
detailed calculations describing this effect for the case of 
two inclusions, and have presented how to approach the 
problem for additional inclusions. 

Wc have considered here only the case of zero total 
topological charge brought about by parallel boundary 
conditions at infinity. If the c-director has some fixed 
angle relative to the outer boundaries of the sample, then 
the total topological charge is +f . In this case, a single 
inclusion will not force the nuclcation of a companion de- 
fect. Each additional inclusion will however, and there 
will be (n — 1) defects in a sample with n inclusions. 
This is entirely analogous to the situation in multiple ne- 
matic emulsions where water droplets are captured inside 
spherical nematic drops 
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APPENDIX A: AN ELECTROMAGNETIC 
ANALOGY 

We take the most general $ to be of the form 

$ = ^ (gf log |r - r,| + qfct>r, (r)) . (Al) 

i 

It can be shown that 

J (V log |r - r^l) • (V</.r, ir))(fx = (A2) 

and thus if we define 

$^ = ^5flog|r-r,| $^=^gf</>r,(r) (A3) 

i i 

then we have 

J {V^f(fx = J {{V^^f + (V$^)2) d^x (A4) 

and we see that the two families of singularities do not 
interact. 

Now recall that the potential for a line charge parallel 
to the ^— axis is 

Alog(V(x-xo)2 + (y-yo)2) (A5) 

where (xo,j/o) are the {x,y) coordinates where the line 
pierces the xy— plane. Thus, we can identify with a 
collection of line charges or equivalently 2D electrostatic 
charges. Next, noting that for any function f{x,y) we 
have 

V/ • V/ = (V X /e,) • (V X fe,) (A6) 

we can define A = $^6^ and note that 

(V X A)^ = (V$^)^ (A7) 

Interestingly, the vector potential for a line of magnetic 
charges parallel to the 2— axis piercing the xy— plane at 
(a;o,2/o) is simply 

A = tan-i ( e, (A8) 

\X-XqJ 

which corresponds precisely to the singularities in 
Equivalently we can view this as a magnetic charge in 
2D, just as we viewed the electrostatic line charge as 
a charge in 2D electrostatics. Now, making the final 
identification, 

E = -V$^ B = V X A (A9) 

we can write 



F='j- / (e2 + b2) d^x (AlO) 
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which looks just Hke the energy density for electromag- [15] Pindak, R., and Moncton, D., 1982, Phys. Today, 35, 5, 
netism. Furthermore, we can define electric and magnetic 57. 
charge densities 

p"" ^Y.if^^^-'^^) p"" ^Y.'ii^^^-^i) (All) 

i i 

and then complete the electromagnetic analogy by writ- 
ing 

E = -V$^ B = V X A (A12) 
V-E = 27rp^ V-B = 27rp^. (A13) 

This is why we have chosen to call the logarithmic sin- 
gularities "electric charges" and the tan^-'^ singularities 
"magnetic charges" . Of course, in this analogy we have 
no reason to expect the lines of magnetic charge to be 
quantized, whereas in the original problem there is a 
clear explanation as to why the tan~^ singularities are 
quantized (the requirement that the physical field $ be 
single- valued). 

Alternatively we could have exploited the fact that 

(V0r„(r))2 = (Vlog|r-ro|)^ (A14) 

treated the tan~^ singularities as a separate species of 
logarithmic singularities and viewed our free energy as 
being the sum of two distinct non-interacting electro- 
static energies. 
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Caption Fig.l 

(a) A schematic overhead view of a generic circular inclu- 
sion, showing the relative positions of the physical and 
virtual defects, as well as the relationships between the 
far-field director, the director at the boundary of the in- 
clusion and the direction of the topological dipole mo- 
ment. The director c is shown at the boundary and its 
angle with respect to the tangent vector t measured in 
a counterclockwise fashion is (j)a- The direction of the 
dipole moment p is given by the usual convention (i.e. 
from negative to positive charges) . We see from this and 
figure (2) that p makes an angle oi — (j)a with re- 
spect to the far-field director, (b) An oblique view of a 
schematic representation of a film containing an inclu- 
sion consisting of one additional layer. The rods repre- 
sent the actual nematogens in the film which are tilted 
with respect to the normal to the smectic planes. On 
the top plane the projection of the nematogens into the 
smectic plane is represented using the conventional nails 
with heads representation, where the head of the nail 
distinguishes which end of the projection into the plane 
corresponds to the end touching the plane. The com- 
panion defect is centered at the dark circle. The figure 
to the right indicates the two dimensional interpretation, 
consisting of the boundary of the circular inclusion, the 
direction of the far-field director as well as the orienta- 
tion of the director at the boundary of the inclusion, the 
positions of the physical and virtual defects and the re- 
sulting topological dipole moment, (c) Photograph of a 
four layer freely suspended smectic-C liquid crystal film 
with a circular five layer island of radius 75 microns. The 
companion —1 defect is seen sitting outside the island. 
The collection of image charges is shown in the inset. 
The presence of only a single -1-1 point defect in the cen- 
ter of the island in the photograph reminds us that the 
charges that the companion —1 defect interacts with are 
image charges imposed by the boundary conditions and 
not real physical charges inside the island. As mentioned 
in the paper, the rigid boundary conditions along the 
dislocation separating the bulk from the interior of the 
island allows us to consider the regions independently. 
The visible — 1 defect has been observed undergoing large 
fiuctuations in its radial position as well as its azimuthal 
orientation. 

Caption Fig. 2 

Three examples of the director field c around an inclu- 
sion for three different values of (j)a- The direction of p 
is indicated. In each case we have placed the external 
defect at its preferred position, in particular Scjjd — and 
hence 4id = 4>^d- values of (f)d, the azimuthal position 
of the companion defect, are indicated. Fig. |^ shows 
examples of fields where 5(j)d is nonzero. Also note that 
in each picture the far-field director is taken to be in the 
+x direction. 



Caption Fig. 3 

Director configurations with nonzero values of 5(j)d- The 
far-field director is assumed to be in the positive x direc- 
tion, and (/)a = 0. Thus the equilibrium position of the 
companion defect is above the inclusion (0^ = 7r/2) and 
the equilibrium direction of p is down. In (a) 5(j)d — 71/4, 
and in (b) 5(j)d = tt. Compared to the equilibrium con- 
figuration of Fig. 2(a), note the increasing distortion of 
the c-director as 54)d increases. The energy of (b) with 
p pointing up is clearly higher than that with p down. 
Unlike 3D nematics, there is a preferred direction rather 
than a preferred axis of orientation for p. This is not 
surprising since a nematic is invariant under n ^ — n 
whereas a smectic-C is not invariant under c — s- — c. 

Caption Fig. 4 

Above we have two interacting inclusions. The equilib- 
rium configuration for the dipoles is represented by the 
solid arrows. The dashed arrows represent possible orien- 
tations of the dipoles due to thermal fluctuations. Note 
the correlations in the directions of the dipole moments, 
when pi is tilted downward we expect to find p2 tilted 
upward. 
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FIG. 1. 
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FIG. 4. 



14 



